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Abstract 

We define compressive and rarefactive waves and give the difTerential 
equations describing smooth wave steepening for the compressible Euler 
equations with a varying entropy profile and general pressure laws. Using 
these differential equations, we directly generalize P. Lax's singularity 
(shock) formation results in [9] for hyperbolic systems with two variables 
to the 3x3 compressible Euler equations for a polytropic ideal gas. Our 
results are valid globally without restriction on the size of the variation 
of initial data. 
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1 Introduction 

In this paper, we consider the initial value problem for the compressible Euler 
equations in Lagrangian coordinates in one space dimension, 

rt-u, = 0, (1.1) 
ut+p, = 0, (1.2) 

(1^2 + e)t + (up), = 0, (1.3) 

where p is the density, t = is the specific volume, p is the pressure, u is the 
velocity, and e is the internal energy. For solutions, (11.31) can be equivalently 
replaced by 

St - 0, (1.4) 



* chen@math.psu.edu 

tThis paper is part of the first author's Ph.D dissertation at the University of Mas- 
sachusetts, Amherst. 



1 



where S is the entropy. When the entropy is a constant (isentropic fluid), (ll.ip 
and (11.21) become a complete system, known as the p-system [15] . 

The formation of shock waves from smooth initial data is one of the central 
problems of conservation laws. The formation of shock waves takes place when 
one or more gradients blows up. In this paper, we consider the formation of 
such kind of singularity from smooth initial data. In [5], Lax gave singularity 
formation results for hyperbolic conservation laws with two variables, including 
the p-system as an example. This says that the singularity forms from smooth 
initial data if and only if the initial data includes compressive waves. For the 
conservation laws with a coordinate system of Riemann invariants, similar sin- 
gularity formation results hold, c.f. [1]. However, for the conservation laws 
with more than two variables and without a coordinate system of Riemann in- 
variants, including (|l.ip ^ (|1.3p . the singularity formation results are only for 
restrictive initial data, see {7^, |T2], |TT]. In this paper, we directly general- 
ize Lax's singularity formation results for the p-system to (|l.ip ^ (|1.3p . The 
singularity formations for the compressible Euler equations in multiple space 
dimensions are considered by [TJ and |T3] . 

Comparing to the p-system, one of the main difficulties for studying (11.11) ~ 
(jl.Sp . even for the smooth solutions, is the smoothly varying entropy field. Wave 
propagation in a smoothly varying entropy field is still not well understood: even 
some basic questions, like what should be considered as the compressive and 
rarefactive waves have not previously been answered. Furthermore, the study 
of smoothly varying entropy field is also the basis for studying the formation of 
singularity, and interactions between shock waves and simple waves. 

In an isentropic domain, people have long known how to define the com- 
pressive and rarefactive waves for smooth solutions by the changes of density, 
pressure, wave speed or Riemann invariants, c.f. [3], [5], [S]: [S], [IS]- How- 
ever, in a smoothly varying entropy field, there are no Riemann invariants, and 
the forward and backward waves can't be geometrically divided. In order to 
define rarefactive and compressive waves, we need to find a variable which is 
not impacted by the variation of entropy but still can discriminate rarefactive 
and compressive waves. Pressure is an appropriate variable, since it is invariant 
inside the stationary entropy field. 

We use superscripts / and \ to denote the directional derivatives on the 
forward and backward characteristic lines, respectively: 

/ = 9( + cdx, \ = 9t - cdx, (1.5) 

where c is the wave speed, c.f. [9]. We define the Rarefactive (R) and Compres- 
sive (C) waves (or R/C character) as follows: 

Definition 1.1. Consider a solution of (|l.ip ^ (|1.3p . which is smooth in an open 
set U in the (i, x)-plane and A is a point in U. We say that the solution is 
forward (backward) rarefactive at A, if and only if p^ < (p' < 0); it is forward 
(backward) compressive at A, if and only if p^ > (p' > 0). 

Here, we use the derivative of pressure in the opposite direction to define the 
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R/C character, which helps us discount the disturbance from the waves in the 
opposite direction. This definition gives us appropriate physical explanation of 
R/C characters. This definition is not only for the simple waves [T5], but also 
for waves in the wave interaction regions. 

Then we restrict our consideration on the polytropic ideal gas dynamics, 
where 

p^Ke^T-^, (1.6) 

K and are positive constants, and 7 > 1 is the adiabatic gas constant. We 
introduce variables z and m: 

z = — -T 2 ^ rn = e^"^ . (1.7) 

7-1 

In this paper, the smooth solution means that u, r, S are C^, which is equiv- 
alent to that u, z, m are and z, m are both positive, by (|1.7|) . If the initial 
data are smooth, m is smooth and positive until either the r or m profile has 
some singularities (gradient blowups) , since the entropy is stationary in solu- 
tions by (|1.4p . Since we only consider smooth wave propagation and singularity 
formation from smooth initial data in this paper, the entropy profile is always 
smooth and stationary. 

We introduce variables 

a = = Ua:+mZ:c + ^mxZ, (1.8) 

13 = = — mzx — ^^-^m^z. (1-9) 

a and /3 are the generalization of Sx and rx in a smoothly varying entropy field, 
where s, r are Riemann invariants in a constant entropy field. 

Theorem 1.2. The smooth solutions in M.l]) ^ i f77g|j satisfy 

a' = A:i{/c2(3a-|-/3)-|-a/3-a2}, (1.10) 

and 

^ki{-k2{a + i(3) + al3~ P^}, (1.11) 

where 

ki — ^2(^-1)" '^"'^^ ' ^2 = 'y(-/+i) -^"^z; "is a positive constant. (1-12) 

The equations (|1.10p and (|l.lip arc not pure ODEs because they aren't closed. 

Then we transform (|1.10p and (|l.lll) into "decoupled ODEs" by changing a 
and /? into new variables y and q, where 

y^ni 2(3t-i) 2;2(-,-i) ((y ^2)2, - 3::pj-ma;z), (1-1'j) 
q = m 2(3t-i) z2(7-i) ((tt - mz):r + j^'Tij^z). (1-14) 
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Similarly, we define 

7 + 1 

y = zW7^({u + mz)^- ^;2^m^z), (1.15) 

q = z5T^ ((u - mz)^ + -^^m^z). (1.16) 

Theorem 1.3. The smooth solutions in i fTTg)] satisfy 

y' ^ ao + a2y^, (1.17) 

q' = ao + a2g^ (1.18) 

where 

02- -"C2(:;^ 



/^e ^-fe^ r 7^-1^ (37+l)(7-l) ^21 

Kc J_i) m^(^^^z^(^~^ < 0. (1.20) 



Clearly, the functions ao and 02 only depend on the density and initial entropy 
profile. These "ODEs" generalize Lax's "ODEs" for p-system in 9 . In fact, 
when entropy is constant, ao = 0, and we recover Lax's equations. In this paper, 
"— denotes the limit when (t,x) approaches a fixed (finite) point {to,xo)- 

For later reference, we state two assumptions that will be used in some of 
our results. 

Assumption 1 : z is not equal to zero or infinity for any x and t. 
Assumption 2 : There exist positive constants Zl, Zjj, AIi ^ M4, such that 
Zl < z < Zu, for all x and t; (1.21) 

Ml < m° < M2, |m°| < Ms, \m°J < M4, (1.22) 

where m^{x) = m{0, x) is given by the prescribed entropy profile s(0, x) = s(t, x) 
by and (T^. 

By ()1.7p , Assumption 1 means that the density is not zero or infinity, and (|1.2ip 
means that the density has positive upper and lower bounds. Furthermore, all 
realistic smooth (C^) initial entropy profiles satisfy (|1.22p . We impose Zl to 
avoid potential problems at vacuum, which are addressed in an upcoming paper 
[20] . Zl can be arbitrarily small. 

Using (|1.17l) and (|1.18p . we can now give some singularity formation results. 

Theorem 1.4. Assume the initial data are smooth, and Assumption 2 holds. 
The constants Zjj and Mi^M^ are given in Assumption 2. Then there exist 
positive constants N and N depending only on Zu and Mi^M^, such that, if y 
or q is less than —N or y or q is less than —N somewhere in the initial data, 
then \ux\ and/or \tx\ blow up in finite time. When 7^3, all the results hold 
without the lower bound of density in Assumption 2. (Here y, q, y and q are 
defined in f7jg)) ^ fr7g|) . ) 
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Theorem 1.5. Assume the initial data are smooth, and Assumption 2 holds. 
Furthermore, assume that there exists a point A G 'M., such that the initial 
entropy profile satisfies 

^ 0, for x>A. (1.23) 

If, in the initial data, yo — y{0,A*) < at some point x = A* with A* > A, 
then \ux\ and/or \tx\ blow up before some finite time T*, where 

0<T, sc: -, (1.24) 

yo min(-a2) 

and min(— 02) is a positive constant depending on {Zjj) o,nd Mi {M2) when 
1 < 7 < 3 (7 > 3). When 7 3, all the results hold without Assumption 2. 
Symmetric results hold for q. (Here 02 is defined in il.20\) .) 

Theorem 11.41 implies that gradients of solutions blow up if the initial com- 
pressions are strong enough, which is a direct generalization of Lax's singularity 
formation resuhs for p-system. When the variation of entropy is mild, N and 
N are close to zero, so the shock free solutions are "almost rarefactive" , which 
is consistent with Lax's singularity formation results. In |16) and a forthcoming 
paper [5], examples are given, of solutions containing compressive waves, but 
the gradients of the solutions don't blow up. 

Our "ODEs" and the singularity formation results are for arbitrarily large 
smooth initial data, where there is no restriction on the amplitude of the waves. 

We divide this paper into 7 sections. In sections 2 and 3, we review the 
background and definition of R/C character in the p-system. In section 4, we 
define rarefactive and compressive waves in a smoothly varying entropy field. In 
sections 5 and 6, we give "ODEs" for smooth solutions. In section 7, we prove 
the singularity formation results. 



2 Equations and Coordinates 

In this paper, we focus on the polytropic ideal gas dynamics, where the equation 
of state is given by, 

PT = RT, (2.1) 

and 

e = c„T = ^-pr, (2.2) 

7- 1 

with 

p^Ke^T-^. (2.3) 

Here S is the entropy, T is the temperature, R, K, Cy are all positive constants, 
and the adiabatic gas constant 7 > 1, c.f. [3]. The (Lagrangian) wave (sound) 
speed is given by 

c = V-Pt = A/^r"T~e~. (2.4) 
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We use the coordinates provided by B. Temple and R. Young, in [TC]. Define 
new variables m and z for S and r, by 

s 



and 



(2.5) 



-dT=^— ir- — , (2.6) 
where we use (|2.4I) and (|2.5p . It follows that 

T = if^2~^, (2.7) 

p = Kpm'^zT^, (2.8) 

c = c{z,m) = Kcinz-i-'^ , (2-9) 
where i^T, -?''p and Kc are positive constants given by 

Kr = (^)^, (2.10) 

7-1 

Kp = KK-\ (2.11) 

Kc = y^X;"^, (2.12) 

and 

Kp - ^^i^c. (2.13) 

27 

By (|2.5p and (|2.6I) . for solutions, the Lagrangian equations (|l.ip ^ (|1.3p are 
equivalent to 

zt + —u^ = 0, (2.14) 
777 

Ut + 777.CZ2, + 2 — 777^ = 0, (2-15) 
777, 

777t = 0, (2.16) 

where the last equation is coming from (|1.4p instead of (|1.3I) . c.f. [15| . 

When the entropy is constant, p — p{t), so (|2.14p ^ (|2.16p change to the 
p-system: 

Zt + —u, = 0, (2.17) 

777 

Ut + 777CZ2; = 0. (2-18) 

So the corresponding Riemann invariants are 

r = u — mz, (2-19) 

s^u + mz, (2.20) 

which satisfy, by (PTT)) and ([^T^ . 

n - cr, = 0, (2.21) 

St + cs^ = 0. (2.22) 
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3 Rarefactive and Compressive Waves in a Con- 
stant Entropy Domain 

We first review the definition of tlie rarefactive and compressive waves for 
smootli solutions in a constant entropy domain, where the system is 2 x 2. At 
this time, the system has full set of Riemann coordinates (|2.19p and (|2.20p . For 
a single simple wave, the wave is rarefactive (compressive) if pressure, density 
or wave speed, are decreasing (increasing) from ahead of to behind the wave, c.f. 
[8]. Here we consider not only the simple waves, but also the local rarefactive 
and compressive characters at some points or open sets in the wave interaction 
regions, c.f. [9], [16], [19]. Sometimes, we use the names "rarefactive wave" 
and "compressive waves" or even "i?" and " C" , instead of the rarefactive and 
compressive character at some points or open sets, without confusion. 

In this section, we assume entropy is a constant, hence p is a function which 
only depends on r, so that 

p'=PrT\ c'^CrT\ p' = -\t\ (3.1) 

Recall subscripts "/" and "\" denote the directional derivatives along forward 
and backward characteristics, respectively, which are defined in (|1.5p . Note 
Pr < 0, and Ct- < since Prr > and (|2.4p . Furthermore, 



£t 

c 



(3.2) 



where we use (P^ . ^^TM and pTT]) . respectively. 
Hence 



t' < ^ c' > ^ p' > ^ s,; < ^ St > ^ p' > 0, (3.3) 

which means that the forward waves are compressive; and 

t' > ^ c' < p' < > ^ St < ^ p' < 0, (3.4) 

which means that the forward waves are rarefactive. 

In fact, (|3.3p and (|3.4p mean that pressure, density and wave speed, are 
locally increasing or decreasing from ahead of to behind the forward wave along 
backward characteristic line, respectively. We use the directional derivative 
along backward characteristic to define the forward R/C character, which dis- 
counts the disturbance from the backward waves. Symmetrically, we define the 
backward R/C character. 



4 Rarefactive and Compressive Solutions in a 
Smoothly Varying Entropy Domain 

In the previous Section, we describe the equivalent definitions of rarefactive and 
compressive waves in a constant entropy domain. However, when the entropy is 
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smoothly varying, there are no Riemann invariants anymore, and the forward 
and backward waves can't be geometrically divided, which give difficulties for 
defining the R/C characters. In this Section, we present a definition of rarefac- 
tive and compressive waves in a smoothly varying entropy domain, which is new 
and fundamental for the smooth waves of the compressible Euler equations. 

The conditions shown in the previous Section are not equivalent anymore 
in the domain with smoothly varying entropy profile. This is because pressure 
and wave speed become functions depending on both density and entropy. In 
order to pick up the right condition, we first consider the stationary solutions of 
(|l.ip ^ (|1.3p . where u, t, S are stationary. So u, p are constant by (|l.ip ^ (|1.3p . In 
particular, the density need not be constant when the entropy profile is varying. 

We go back to check the rarefactive and compressive conditions given in 
the previous Section. In the stationary solutions, p' = — while the other 
quantities are nonzero if the entropy profile is smoothly varying. In the other 
words, pressure is the only thermodynamic variable which is not impacted by 
the stationary entropy field. Moreover, change of pressure contributes to the 
wave propagation on the (non-vertical) characteristic lines, i.e. rarefactive or 
compressive waves. So pressure is the only possible variable which still can 
distinguish rarefactive and compressive waves in a smoothly varying entropy 
field. So we can get the definition of rarefactive and compressive waves in 
Definition 11.11 Note that we use the directional derivative of pressure along the 
opposite characteristic to define the R/C character, which helps us discount the 
disturbance from waves in the opposite direction. This definition also applies 
for the general pressure law with pr < 0, prr > 0. 

We define 

P' . P' 
a = — -, /3 = — -. 
c c 

Lemma 4.1. 

a = Ux + mzx + ^^rrixZ, 

and 

j3 — Ux ~ mZx — ^^^^rrij-z. 



(4.1) 

(4.2) 
(4.3) 



Proof. 



-c a = p 



{K^virzi-^ y 



KpTn'^^^z-'''^ zt — cKpTTi? -^^^z-i-^ Zx — 2cKpmmxZ-<-^ (4.4) 



-Kpmc^^z~ {ux + mZx + ^nixz) 



= c^{ux+mzx + ^mxz), 



where we use (|2.9p and (|2.13p . Similarly, we can prove (j4.3|) . □ 

By the Definition ll.il and (|4.1I) . we can equivalently define the R/C character 
by a and /3. 
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Lemma 4.2. In the polytropic ideal gas, the local R/C character of the smooth 
solution is given by: 



Forward 
Forward 
Backward 
Backward 



When Assumption 1 holds, 



R 


iff 


a > 0, 


C 


iff 


a < 0, 


R 


iff 


/3>0, 


C 


iff 


/3 < 0. 



oo iff 



or oo. 



(4.5) 



(4.6) 



Proof. Clearly 
and, 

By g3 and (li3D . 



a + /3 ~ '2ux, 
a — j3 — 2(mZx 



7-1 



rrixz). 



(4.7) 
(4.8) 

(4.9) 
(4.10) 



By (112]), (1131), dlH), dmni), (EH) and Assumption 1, we get (gSl), where we also 
use that m is stationary and positive in smooth solution by (|2.5p and (I2.16p . □ 

Note: We give the definition of compressive and rarefactive waves in a 
smoothly varying entropy profile from a purely physical point of view, which also 
explains the physical meanings of Riemann invariants in the p-system. In fact, 
in the p-system, the physical meanings of the derivatives of Riemann invariants 
can be explained by the directional derivatives of pressure: 



P 



P 



(4-11) 

where we use that m is a constant, (I2.19p . (j2.20p . (|4.ip and Lemma |44] So a, 
/3 can be considered as the generalization of Sx and r^, in a smoothly varying 
entropy field. When = 0, a and /3 equal to Sx and rx, respectively. 



5 Smooth Wave Propagation 

In this section, we consider the smooth wave propagation of (11.11) ^ (11.3^ . By 
considering the directional derivatives of a and /3, we construct "ordinary dif- 
ferential equations" (|1.10p and (|l.lip in Theorem 11.21 which together with the 
definition of rarefactive and compressive waves will give us a framework for 
the wave propagation of smooth solutions. These "ODEs" are not real ODEs 
since the derivatives are in different directions and (jl.lOl) and (jl.lip depend on 
the variable z, but they are much simpler than the original partial differential 
equations. 
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5.1 Proof of Theorem [HI 

Proof. By Lemma l4.ll 

a' = (ux + mzx + '^^mxz)' 

= u^t + mz^t + ^m^zt + c[u^.^ + mz^^ + m^z^ + ^{m^z^ + m^^z)] 
= (u^t + cmz^a:) + (cMxx + mz^t) + ^^cm^z^ + 2^m^Zi + ^cm^a^z. 

(5.1) 

By dm])- dsns]) and 



(wt + cmza;)a; = = -(2i4rprnma;ZT-i (5.2) 
m 

and 

(cMa; + mzt)a; = 0. (5.3) 

So 

Ua:t + cmZxx = —2KrjmnxZi~^ Zx — ^^^^Kc'm'^z-i-'^ i^x)'^ 

-2Kp{mx) z-i-^ - 2KpmmxxZ-'-^ - ^i^Kp-mnixZ-i-^ Zx, 

(5.4) 

and 

CUxx + WZa;t = -CxUx - TUxZt = -^^KcmZ~ ZxUx- (5.5) 

Furthermore, the sum of the second term in the right hand side of (|5.4p and the 
right hand side of (15.51) is 



-^Kcmz^-i Zx{mzx + Ux) 
= —^^^^KcZ^i~^ {a — ^^^rUxZ — Ux){ce — -^^Tnxz) 
= -^KcnixZ-'-'^Ux - ^-::^Kc{mxyz-<-'' 

+ X±l^K,z^ [{^fuxz + ux)a - a% 

By and ((^ . the right hand side of ([5TT|) equals to, 



(5.6) 



'2Kcmmx z-<-^ Zx — 2Kp{mx) zi-'^ — 2KpmmxxZ''-^ — ^^^KpUirnxZi-'^ z. 

— — — KcXnmxZ-i-^ Zx '-^KcfTixZ-i-'^Ux + -^z^KcmnixxZ''-^ 

+ ^Kr^z^[{'^mxZ + Ux)a ~ a 

(5.7) 

where we use (|2.14p to get rid of zt- Plug (|2.13p into (|5.7p . then (15.71) equals to, 
-2ii:c'Ti;jZ~U3; - 27 -'y-i j^^(m^)2;^7^ _|_ zi^Xzl ^ ^mm X z^ z : 



-1+1 

+ 2±l^K,z-i^^[C-^mxZ + Ux)a-a\ 



-y — 1 L\ -y 

Using mzx = a — Ux — -^—^rnxZ, which is from (|4.2p . (|5.8p can be simplified to 
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i^cz^{im.^i. + (im, + ^^)a-^^}. (5.9) 
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By dm, 

Ux = (5.10) 

Plugging ([CTT)) into ([O)) . we get ([TTTU|) . The calculation of PTTT]) is same. □ 

Corollary 1. In an open set U in the (t, x)-plane, assume the solutions of 
lll.l]) ^ JTTSj) are smooth, nrix 7^ 0, and Assumption 1 holds. If /3 = in U , then 
a = in U ; if a = in U , then /3 = in U . 

Proof. If /3 = in ?7, —kik2a = by (jl.lip . hence a = by Assumption 1. 
Another case can be proved similarly. □ 

In this corollary, we show that there will be no "pure" forward or backward 
rarefactive or compressive waves inside any open set in (t, a;)-plane with varying 
entropy, i.e. the forward and backward waves will exist or die out mutually. 
This is different from the isentropic domain, where we can find rarefactive or 
compressive simple waves. 

6 "Decoupled" ODEs 

Next, after introducing new variables y and q, we change (|1.10p and (II. lip into 
"decoupled ODEs" (fTTTl) and ([TTTSl) in Theorem [TH which are only coupled 
by z and x. The variables y and q are given in (|1.13p and ()1.14p . These new 
"ODEs" will help us prove the singularity results. The equations (|1.17p and 
(jl.lSp generalize Lax's "ODEs" for p-system in [5]. In fact, when the entropy is 
constant, oq = 0, so (|1.17p and (|1.18p change to 

y' = a2y^, q' = a2g^, (6.1) 

which is exactly the "ODEs" provided in ^ for the p-system. 

6.1 Proof of Theorem [HI 

Proof By (l2T4l) . (EH) and (l43l) . 

z' = Zt+ CZx 

= - — {ux-mzx) 
m 

= -K,z-^^{P+^^mxz). (6.2) 
7 



Hence 



P = ^rqrr " - — ^mxZ. (6.3) 

Krz— 7 



Plugging (1131) into ([LT0|) . we get 



a' = fci{fc2(3a 4t:- - ^mxz) + a{ 4n- - ^ruxz) - a^}. 



(6.4) 
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We move the terms including z to the left hand side, then we multiply by 

-y + l 

22(^-1) on both sides. After simplification, we have 

T+l + 1 I , T + 1 , 

Z2(T-1) a' + ^1 _|_ _^tti_Q,2;2(T-l) ^z' 

. 3(7 + 1) 1 3(7 + 1) ^, , 3(7 + 1) -I „ 

(6.5) 

where we use (|1.12p . The left hand side of (|6.5p is equal to 

(a^^ + _^^5(?7^+i^^)'__^^5(?^+i(^^)'. (6.6) 
We define a new variable y, 

T+i 7—1 T+l I 1 

y = q;z2(7-i) -I -! -z2(-'-i)"^ m^,. (6.7) 

7(37 - 1) 

So 

7+1 7 ~ 1 

a = yz 2(7-1) — —zm^. (6.8) 

7(37-1) 

Hence, by (|6.6p . (|6.8p and (m^j)' = cnixx, (|6.5p changes to 

3(7+1) , 1 3(7+1) , , 

9--,, 3(7+1) 7 + 1 

+ (yz _ _2_^2;^^) 

-j^Kj^-\yz~wrTj -^^^zm^y. (g.g) 
After simplification, we have 

y' ^ do + aiy + a2y'^ , (6.10) 

where 

1,7-1 (37+l)(7-l) 3(7+1) I 1 

ao = K.-[^mm^^- ^ ^ ml]z^^\ (6.11) 

7+1 "> + ! 1 

02 - -KcTTT -Z^i^-\ (6.13) 

2(7 - 1) 

Then we do one more simplification by multiplying 

3(3-7) 

l_i = m 2(37-1), (6.14) 
on ()6.10p . In fact, it is easy to check that 

jl' = —5,1/2, (6.15) 
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y = My- (6-16) 

Hence (|6.10p changes to 

y' = ao + a2y^, (6.17) 

where 

ao = p,ao, a2 = a2/M- (6.18) 
Similarly, we prove prT5| . By (PHI) . (I2H) and 



= Zt ~ CZx 



c 



m 



Hence 



7 + 1 -7 — 1 

-KcZ~{a-- m^z). (6.19) 

7 



7-1 

^,+1 I m^z. (6.20) 



^, 1 3(t+1) I -I „ 9 ^, 3(7 + 1) . 3(-f + l) , 



Kcz— 7 
Plugging (lOO]) into ([m]) . we get 

13' = fci{-fc2(3/3 + ^m,z) + /?( + 2:-im,.z) - p^}. 

We move the terms including to the left hand side, then we multiply by 

7 + 1 

^2(7-1) both sides. After simplification, we have 

0,0-1 3(7 + 1) -I 

(6.22) 

where we use (|1.12p . The left hand side of (|6.22l) is equal to 

{Pz^^ - -^^z^+^rn^y + _LzL_zWr^)+\m^)\ (6.23) 
We define a new variable g, 

7 + 1 7—1 7+1 I 1 

q = Pz^^^ -z^T^+m^. (6.24) 

7(37- 1) 

So 

7+1 7—1 , ^ 

^ = + — i -zm^. 6.25 

7(37- 1) 

Hence (|6.22p changes to 

o, 1 3(7 + 1) I 1 3(7 + 1) I -1 

9 = ^(;^.y_i) j^c^^'-''^'"^ mm^^ - X-^X^z; 2(7-1) 

9_ 3(7+1) 7 + 1 _. 

-^K^m^z^i-r-^){qz ^c-^) + 7. ' zm^) 



27 ^^c'-x^ ^y-^ -I j(3j-i) 

-1 3(7+1) -1 7+1 -1 

2f^Kc^^-'(g^"^ + ^(i^zm,)2, (6.26) 
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where we use (rrixy = —cnrixx- After simplification, we have 

= ao - ai9 + 029^, (6.27) 
where hi are in (|6.1ip ^ (|6.13p . Similarly, we denote 

q = M, (6.28) 

where p, is in (|6.14p . Since 

jl' = ai/2, (6.29) 

(|6.27p changes to 

= aQ + a2q^, (6.30) 
where oq, 02 are defined in (|6.18p . □ 
It is easy to check that a2 < 0, 

ao I ^ (87 - l)mm,, | (87 + l)ml ^ (m"^),, = 0. (6.31) 

So the sign of ao only depends on the entropy profile, for fixed 7. Note that the 
entropy profile of smooth solution is stationary because of (|2.16p . So the sign 
of flo only depends on the initial data. 

Corollary 2. When Assumption 1 holds, 

\y\ or \q\ ^ 00 iff \ux\ or \tx \ co. (6.32) 

Proof By (gJl), (ICTH)) . (jOIl and (lOSl) . 

3(3-7) t+i 

y + q = 2m 2(3t-i) u^-z^c^-i) , (6.33) 

3(3-7) , 1 + 1 37—3 7 + 1 I 1 

= 2m '^(■'-'-''Umz'^'--'-'') -\ ■ z2(-'-i)+ m:^)- (6.34) 

37 — 1 



By the same argument as the proof of (j4.6p in Lemma 14.21 (|6.32l) is right. We 
complete the proof. □ 

Similarly, when Assumption 1 holds, 

\y\ or \q\ 00 iff \ux\ or \tx \ 00. (6.35) 



7 Singularity Formation 

In this section, using the quadratic "ODEs" in the Theorcm ll.31 we will give sev- 
eral singularity formation (gradient blowup) results for ()l.ip ^ (|1.3p with smooth 
initial data. By Corollary[2l we need to find conditions on the initial data, under 
which \y\ or |g| blows up in finite time. 

We consider the quadratic equation (^ = y or q), 

0-00 + 02^', (7.1) 
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where uq and 02 are defined in (|1.19p and (|1.20l) . When oq < 0, there is no real 
root; when ao = 0, there is a unique root 0; when qq > 0, the two roots of (|7.1I) 
are 

±V/^ - ±V .(.+TK3?-i) ("^"^--i^"^^) z^'-'m-^. (7.2) 

By considering ()1.17p and (jl.lSp , we summarize the dynamic system properties 
for smooth solutions in Fig. [1] Moreover, {y > 0, g > 0} is an invariant domain 
for the region with initial entropy profile satisfying oq ^ 0, by Fig. [T] 




Figure 1: ao > 0, oq = 0, ao < from left to right. The arrows indicate y or q 
increases or decreases. 

For the p-system, m is a constant. Hence, ao = 0, and 

3(3-1) 1+1 3(3-1) 1 + 1 

y = Sxm 2(3t-1)^2(^7^^ q^^^^ 2(3^-1)^2^7^^ (7.3) 

so y and Sx [q and r^) have same signs. Hence, the p-system follows the middle 
picture of Fig. [TJ When y or g is negative somewhere, y or q will approach 
negative infinity in finite time. When y and q are nonnegative everywhere in 
the initial data, there is no such kind of singularity, since > 0, 5 > 0} is an 
invariant domain. These results are given by Lax, c.f. [9|. 

We don't expect such "clean" results when the entropy profile is varying, 
because the dynamic properties in this case are much more complicated than 
those in a constant entropy domain, as we can see from Fig. [1] In [16J and 
a forthcoming paper [2], there are shock free examples with compressive waves 
and entropy jumps. We also expect the smooth examples including compressive 
waves. 

Using the nonlinear Riccati type equations (|1.17p and (jl.lSp . we can prove 
Theorems 11.41 and 1 1 . 51 for breakdown of y and q, which directly generalize Lax's 
results in |^9j . This kind of formation of singularity corresponds to the formation 
of a shock. 

Theorem 11.41 savs that the shock free solutions don't include "strong" com- 
pressive waves, since y, g, y, q and a, (3 are larger than some negative constants 
in shock free solutions. When M3, M4 are close to zero, i.e. the variation of 
entropy profile is mild, N and N are close to zero, where N and N are given at 
l|7.5p and (|7.10p . The shock free solutions should be "almost rarefactive" . The- 
orem [L4l is quite sharp at this time, since it is consistent with Lax's singularity 
results in 1^. In fact, when entropy is a constant, N = N — 0, Theorem 11.41 is 
the same as the singularity results for the p-system in [9]. The bounds TV and 
N only depend on Zu , Mi, M2, M3, M4. N and N don't depend on Zl, so Zl 
can be arbitrarily small. 
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In Theorem ll.5[ we consider the singularity formation for some special initial 
entropy profiles. Here, we give an exact example with entropy profile satisfying 
condition (|1.23p and bounded away from infinity and negative infinity. For 
example, 

satisfies (jl.23p . where f{x) is a smooth function that makes m{x) smooth and 
bounded away from zero and infinity. A large class of entropy functions satisfies 

7.1 Proof of Theorem 11.41 

Proof. By (|7.2p and Assumption 2, the roots of (|7.ip . if these exist, have uniform 
lower bound — iV, where if 1 < 7 ^ 3, 



3(3-7) 



^=(l + e)V .(,+T)Xi) ^W ^^'^-^'^ A/r^'^^-^ (7.5) 

if 7 > 3, Ml changes to M2 in (|7.5p . Here, e is a fixed positive constant, which 
can be arbitrarily small. So 

"0 + "2 ^1 ^ g)2 < 0' (7-6) 

since a2 < 0. Recall, m is stationary by (|2.16l) . 

If J/ < — somewhere in the initial data, then by (|7.6p and 02 < 0, 

«o + a2^^|^<0. (7.7) 



So ([TTTI) gives that 



2/' = 00 + 022/' < (1 - 7— ^)a2y'. (7.i 

(1 + e)2 



Hence y is decreasing at this time by (|1.20p . so y is always less than — A^ along 
forward characteristic line. Furthermore, by (|7.l 



lit'. 1 



»«)^i(0)-/„'^-(TT7F)-'"- <") 

where the integral is along the characteristic. By (|1.20p and Assumption 2, 02 
is negative and bounded below. So the right hand side of (17. 9p approaches zero 
in finite time, which gives that yit) approaches —00 in finite time. By ()6.32p . 
\tx\ and/or \ux\ blow up. 

We have another version of the singularity formation results if we start from 
(|6.10p and (|6.27p for y and q. We assume that 



iV = (1 + ^) (7-l)[ \,-^,\M^\M\M^ + mM. ^ Zu^^\ (7.10) 
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where e is a fixed positive constant which can be arbitrarily small, 



A2 = 1^^^(2472 + 167-8). 



(7.11) 



The solutions of 

5o ± + Y^a2r = 0, (7.12) 

if existing, are 

y^^^l 2(37-l)(7+l) 

which are always larger than —N. If y < —N somewhere in the initial data 

y' < (1 - (r^)^2y'- (7.14) 

Hence y will go to negative infinity in finite time, by the same proof of the 
singularity result for y. By (j6.35p . the absolute values of and/or blow up. 
We have symmetric results for q and q. □ 



7.2 Proof of Theorem fTTSl 

Proof. We use ?/(t) to denote the function y along the forward characteristic 
starting from (0,A*), and assume yg = y(0) < 0. By (|6.3ip . ag ^ when 
x> A. So (|1.17p changes to 

V' > a^y". (7.15) 

By solving (|7.15p . we get 

11/" 

- / a2dt, (7.16) 

y J 

where the integral is along the characteristic. Hence, when the right hand side 
of (|7.16p goes to zero, y goes to —00. So, by i/o < 0, (ll.20p and Assumption 
2, the blowup happens before T, ^ — .^^i^-a^) ■> where min(— a2) is a positive 
constant depending on {Zij) and M\ {Mi) when 1 < 7 < 3 (7 > 3). When 
7 = 3, 02 is a constant, so we don't need Assumption 2. Symmetric results 
apply for g. □ 

When 1 < 7 < 3, we don't need the upper bounds for z, \mx\ and |ma;a;|. 
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